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Abstract. We derive a family of series representations of the multiparame- 

i i ter fractional Brownian motion in the centred ball of radius R in the N- 

0^ dimensional space l". Some known examples of series representations are 

Q, shown to be the members of the family under consideration. 
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1. Introduction 



The fractional Brownian motion with Hurst parameter H G (0, 1) is defined as the 
centred Gaussian process £(£) with the autocorrelation function 

R(s, t) = E£( s K(i) = \{\s\™ + \t\ 2H - \s - t\ 2H ). 

This process was introduced by Kolmogorov (1940) and became a popular statis- 
tical model after the paper by Mandelbrot and van Ness (1968). 

There exist two multiparameter extensions of the fractional Brownian mo- 
tion. Both extensions are centred Gaussian random fields on the space R N . The 
multiparameter fractional Brownian sheet has the autocorrelation function 

1 N 

R (*>y) = ^U(M 2Hi + \Vi\ 2Hi \xi-Vi\ 2Hi ), h j e (o,i), 

while the multiparameter fractional Brownian motion has the autocorrelation func- 
tion 

i?(x,y) = i(||x|r+||y|r-||x-y|r), (1.1) 
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where || • || denote the Euclidean norm in R^ and where H £ (0, 1). 

Malyarenko (2008) derived a series expansion of the multiparameter fractional 
Brownian motion. His expansion converges almost surely (a. s.) in the space C(B) 
of continuous functions in the centred ball 8 = {x£ R N : ||x|| < 1 }. 

In fact, the above mentioned series expansion is a member of a family of 
series expansions. To describe this family, introduce the following notation. Let 



c NH = J2ir( N - 2 )/ 2 T(N/2 + H)T(H + 1) sm(irH), 



(1.2) 



where T denote the gamma function. A review of special functions is given in 
Section 2. Denote 



a m (s,u) 



„,H-i/2r< 



2H-r 



where G 



p. 1 1 



cnhs 



Q,\ , . . . , Q n , G n _|_i , . . . , dp 

b\, ■ . . , b m , b m+ i, . . . , b q 




0. 



m > 1. 



(1.3) 



is the Meijer G-function. 



Let Z + be the set of all nonnegative integers. For a fixed m £ Z+, there exist 

/tKJV)= (^ + iV-2)(m + Ar-3)! 
v ' ; (N-2)\m\ 

different real- valued spherical harmonics 5^j(x/||x||) of degree m. Fix R > 0, and 
for each pair (m, I) with m £ Z + and 1 < I < h(m,N), let {e l mn (u): n > 1 } be 
a basis in the Hilbert space L 2 [0, R]. Let b l mn (s) be the Fourier coefficients of the 
function a m (s,u) with respect to the introduced basis: 



b l mn ( s )= a m (s,u)e l mn (u)du. 



(1.4) 



Finally, let {C mn : m £ 2 1+ ,n > 1,1 < I < h(m,N)} be the set of independent 
standard normal random variables. 



Theorem 1. For any choice of the bases {e mn (u): n > 1}, the multiparameter 
fractional Brownian motion £(x) has the following series expansion 



e(x) 



oo h(m,N) oo 

E E J2 bl mn(M)SL(^/M)Ln- 
m— / — 1 n— 1 



(1.5) 



The series (1.5) converges in mean square in the centred closed ball Br — { x £ 

R N : ||x|| <R}. 

In Section 2 we review the necessary definitions and properties of some special 
functions. In Section 3 we give an outline of proof of Theorem 1. Section 4 contains 
examples, while proofs of technical lemmas are postponed to Section 5. 
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2. Special functions 

This section is intended for readers who are not experts in the theory of spe- 
cial functions. We review definitions and state some elementary properties of the 
relevant special functions. 

This material can be found online at http://functions.wolfram.com, a 
comprehensive online compendium of formulas. 

2.1. The gamma function 

The gamma function of a complex variable z with Re z > is defined by 1 



I» = 
By partial integration, we obtain 2 

T(z 



^e~ l dt. 



1) 



T(z) 



This formula is used to extend T to an analytic function of z € C \ Z_ 
Z_ = {0, —1, . . . , — n, . . .}. The points z e Z_ are the poles. 



(2.1) 
where 



2.2. The Meijer G-function 

Let m, n, p, and q be four nonnegative integers with < m < q and < n < p. Let 
Oi, . . . , a v , b\, . . . , bg be points in the complex plane. Assume that for each k = 1, 
2, . . . , p and for each i = 1, 2, . . . , q we have a^ — 6; + 1 ^ Z + . Then, there exists 
an infinite contour C that separates the poles of T(l — a^ — s) at s = 1 — au + j, 
j € Z + from the poles of T(pi + s) at s = —hi — I, I £ Z + . The Meijer G-function 
is defined by 3 



g: 






,b q 



2ni 



1 Cm, tt n -|_l, . 

! Urn: Vm-\-l i ■ 

T(s + 6Q ■ ■ ■ T(s + 6 m )r(l - tii - s) • • ■ r(l - a n - s) 
T(s + a n+1 ) • • • T(s + a p )r(l - b m+1 - s) ■ ■ ■ T(l - b q - s) ' 



The number p + q is the order of the Meijer G-function. 

The classical Meijer's integral from two G-functions 4 is: 



1 /~1S 



Ci,.. 

di,. 



,Ct,Ct+l, ■ 

, d s , a s+ i, 






v rtm.ii 



Ol, 

6i, 



7 Vjn-i Oyn+l j - 



■ ,a p 

..,bg 



dr 



— a. ^m-\~t ,n+ s 
v+p,u-\-q 



<G 



Ol, 

61,. 



■ : ^mi 



a - 



di, 

Ci,. 



a 
a - 



u v , O n -^-l , . 



. , tip 

■ A 



(2.2) 



: http : //functions . wolfram . com/06 .05.02. 0001 . 01 
2 http : //functions . wolfram . com/06 .05.16. 0004 . 01 
3 http : //functions . wolfram . com/07 . 34 . 02 . 0001 . 01 
4 http : //functions . wolfram . com/07 . 34 . 21 . 0011 . 01 
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A plenty of special functions are specialised values of the Meijer G-function. 
In particular, the Gegenbauer polynomials C^(x) appear as 5 

G o,2 ( a,c \ r(26-2c+2)fl(|z|-l) 

2 ' 2 ^ Z b,b+ 1/2 J T{a-2b+c-l/2)(2(a-2b + c-l)) 2b _ 2c+1 (2.3) 

where (a)„ = T(a + n)/T(a) is the Pochhammer symbol, and 



0(x) 



1, x>0, 
0, x < 



is the unit step function. In particular, the Legendre polynomials P n {x) are just 
Gegenbauer polynomials with A = 0. 

The generalised hypergeometric function v F q {a\, 



q * p[ax > ■••. op. °i. •■•>»«; *J - r(ai) . . . r(ap) S, 9 +i 

The Bessel functions appear as 



a p ; 61, . . . , b q ; z) appears 

1 — ai,...,l — On 



0,1-61, 



l 'p 

,1-6, 



J„(z) — G ; 2 



i//2,-i//2 



(2.4) 



2.3. Spherical harmonics 

Let m be a nonnegative integer, and let mo, mi, . . . , totv-2 be integers satisfying 
the following condition 

m = mo > mi > • • • > niN-2 > 0. 

Let x = (x\,X2, ■ ■ ■ ,Xjv) be a point in the space M. . Let 



Tk 



'k+l 



"fc+2 



L N' 



where k = 0, 1, . . . , N — 2. Consider the following functions 



H(m k ,±,x) = 



(x N - 



ixn 



TN-2 



±m« 



JV-3 



n 



miv-2 TT mk-m k+1 

1 at n I I T t. 



' N-2 11 'k 



xC, 



,n, ,+l.\ -/. -2)/2 /^fc + 1 



m fc -m fc + i 



and denote 

y(m fc ,±,x)=r -" l J ff(m fc ,±,x). 

The functions y(mfe,±,x) are called the (complex-valued) spherical harmonics. 
They are orthogonal in the Hilbert space L 2 (S ,JV_1 ) of the square integrable 



5 http : //functions . wolfram . com/07 . 34 . 03 . 0105 . 01 
6 http : //functions . wolfram . com/07 .31.26. 0004 . 01 
7 http : //functions . wolfram . com/03 .01.26. 0107 . 01 
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functions on the unit sphere S N ~ 1 , and the square of the length of the vector 
F(mfc,±,x) is 

,, , t^ ^ 2 ^ 2mfc - J v+ 2r(mfc _ i + mfc + jy _ i _ fc ) 

fcii (m fc -i + (iV-l-fc)/2)(m fe _i-m fe )![r(m fe + (7V-l-fe)/2)]2- 

Let Z = l(rrik, ±) be the number of the symbol (mo, mi, . . . , mN_ 2 , ±) in the 
lexicographic ordering. The real-valued spherical harmonics, S^x), can be defined 

as 

{F(m fc , +, x)/y/L(m k ), rn N - 2 = 0, 

v / 2Rey(m fc ,+,x)/^i(m fe ), mAr_ 2 > 0, 1 = Z(m fe ,+), 
-\/2Imr(m fe ,-,x)/^i(m fe ), TOjv_ 2 > 0, 1 = l{m k ,-). 

The spherical harmonics S^(x) form a basis in the Hilbert space L 2 (S' Ar_1 ). 

3. Proof of Theorem 1 modulo technical lemmas 

Note that the multiparameter fractional Brownian motion £(x) is weakly isotropic, 
i.e., the autocorrelation function (1.1) is invariant with respect to the group O(N) 
of the orthogonal matrices of order N. Let t > 0, let St = { x £ R^ : ||x|| = t } be 
the centred sphere in the space R , and let dw be the Lebesgue surface measure on 
St- Let 7?(x) be a centred weakly isotropic random field. In fact, the autocorrelation 
function i?(x,y) of the random field r/(x) is a function R(s,t,u) of the three real 
variables s = ||x||, t = ||y||, and u being the cosine of the angle between the vectors 
x and y. Yadrenko (1983) proved that the stochastic processes 

X l m {t)= f r,(x)S l m (x/\\x\\)du 
J s t 

are centred and uncorrelated. The autocorrelation function of the process X m [t) 
is 

fl TO (*,t) = c / fl(s,t,u)C , ^ r - 2 )/ 2 (u)(l-u 2 )( JV - 3 >/ 2 du, 



with 

_ 2 N - 2 Tr( N -V/ 2 mW((N -2)/2) 

C ~ ~ T(m + N-2) ' 

The random field ry(x) can be represented as 

oo h(m,N) 

^( X )=E E ^(l|x||)^(x/||x||), (3.1) 

m=0 1=1 

where the series converges in mean square for any x £ "S. . 

The multiparameter fractional Brownian motion £(x) is Gaussian random 
field. Therefore, the corresponding stochastic processes X l m (t) are Gaussian and 
independent. 
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Lemma 1. The autocorrelation function of the stochastic process X m [t) has the 
form 



Rm(s,t) = < 



4h „2ff r 2,2 I S 
2 * ^4,4 I ^2 



2 b l ^2,2 I t 2 



1-H,0,N/2,1 \ 

, m = 0, 
0, 1 - H, 1 - H - N/2, -H 

< 3.2 

H + l- m, N/2 + H\ 

, m > 1. 
0, 1 - AT/2 -m / 



Recall that a function a(s, u) : (0, oo) x (0, oo) i— > R is called the Volterra 
kernel, if it is locally square integrable, and 

a(s, u)=0 for s < u. (3-3) 

A Volterra process with Volterra kernel a(s, u) is a centred Gaussian stochastic 
process r/(t) with autocorrelation function 

/•min{s,£} 

R(s,t) = / a(s,u)a(t,u) du. 

Jo 

Lemma 2. T/ie stochastic processes X l m (t) are Volterra processes with Volterra 
kernels (1.3). 

Volterra processes are important in the theory of stochastic integration with 
respect to general Gaussian processes, see Decreusefond (2005) and the references 
herein. 

By Lemma 2, the autocorrelation function of the stochastic process X l m {t) 
has the form 

/>min{s,£} 

R m (s,t) = / a m (s,u)a m (t,u)du. 

Jo 

By (3.3), the last display can be rewritten as 

Rm(s,t) = / a m (s,u)a m (t,u)du, s,t€[0,R]. 

Jo 

By definition of a basis in the Hilbert space L 2 [0,R], we obtain 

oo 
R m (s, t) = J2 bl mn(s)b l mn (t), 1 < I < H™, N). 

71=1 

It follows that the stochastic process X l m (t) itself has the form 

oo 

■^■mV') = / j "ran Wsmn) 
n=\ 

where the series converges in mean square. Substituting this formula to (3.1), we 
obtain (1.5). 

We conjecture that (1.5) converges uniformly a.s. 
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4. Examples 

Example 1. Let v be a real number, and let j„.i < j u ^ 2 < ■ ■ ■ < j u .n < ■ ■ ■ be the 
positive zeros of the Bessel function J„(u). For any v > — 1, the Fourier-Bessel 
functions 

¥v.n(u) = ~ r . ^Ju{Ju,nU), n>\ 

■Jv+l\Jv,n) 

form a basis in the space L 2 [0,1] (Watson (1995), Section 18.24). By change of 
variable we conclude that the functions 



e mn\ u ) 



2 If 



-J V (R 1 ju,nU), 



n > 1 



form a basis in the space L 2 [0, R]. 

To calculate b l mn (s), use (1.4). First, consider the case of m = 0: 



Using (3.3), rewrite this formula as 
rV2 



N/2, 1 \ _ V2u 
OA-H ) RJ v+l (J v , n 



J V (R 1 j„ jn u)du. 



,1 ( \ CNHVZ f Hrf D -l' \ r <2. 

b n(s) = -^-j T-. — r / U J V {R Jv,nU)G 2 

n- J v+l\3v,n) Jo 



21 S 2 



N/2,1 
0,1 -H 



du. 



To calculate this integral, use formula 2.24.4.1 from Prudnikov et al (1990) with 
particular values of parameters fc = I = 1: 



u a -'J v {bu)G™£ \ uu 



Ol, . . . , a p 
,b q 



h,- 



du ■■ 



2" n m,r+\ I ^W 



6 Q P+2 ' ? I b 2 



1 — (a + i^)/2, oi, . . . , a p , 1 — (a — ^)/2 
b 1 ,...,b q 



(4.1) 



where a = H + 1, b = R 1 j v , n , m = p = q = 2, r = 0, oj = s 2 , a% = N/2, a 2 = 1, 
b\ = 0, and b 2 = 1 — H . We obtain 



^On(s) = " — EFTT-G. 



T (A \aH +1^4,2 I 2 7 2 



i _ g+l+z; JV -i -I _ g+l-i; 
2 ' 2 ' ' 2 



0,1 -H 

Formula 8.2.2.9 from Prudnikov et al (1990) states 



P:9 



Ol, 






, wq 



Oi, . . . , a p _i 
b 2 ,...,b q 



(4.2) 



If we put v = 1 — H, then (4.2) decreases the order of the Mejer G-function from 
6 to 4. We get 



&(.„(«) 



c NH 2 H + 1 / 2 R H 



ii.i 



4i? 2 



T In \nH+ 1 3' 1 \ S 2 7 2 



0, N/2, 1 
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For further simplification, use the symmetry relation 5 






Ol, 

61 



,...,b q 



•} w p \ rvn,m 



G 



'!■!> 



l-b U ...,l -ba 



1 — <Zl 



1 - a. c 



We get 



0H+l/2r>H I <,2 -2 

"Onl 5 / 1 — T / . \ -H+l "1,3 



h-H{jl-H,n)3l-H,n 

Then, use the argument transformation 9 
with a = 1 to obtain 



a + ai , . . . , a + a y 
a+bi,. . .,a + b q 



AR 2 



1,1 -N/2,0 



y orrm,'.. , ~ 



ai, . . . , a p 
h,...,b q 



CjssrI ' K i> Jl-H,n n l,l I b Jl-H,n 

"1,3 



J 2 - ff (il-tf,„)jf_ + i,„ 4i?2 

and use the formula 10 

iF 2 (ai;6i,6 2 ;z) = f — G{ 3 [ -z 



AR 2 




0,-JV/2,-l 



1 — Ol 
0,1-61,1-62 



with 01 = 1, 61 = N/2 + 1, 6 2 = 2, and z = s 2 j 2 _ Hn /(AR 2 ) to get 

°" l5J J 2 -*(ji-^)jSL 4i?2 r(JV/2 + i) 

x iF 2 (l;iV/2 + l, 2; -s 2 j 2 _ H J{AR 2 )). 
Finally, use the formula 11 



iF 2 (l;2,c;z) = / c _ 2 (2^)r(c)z- c / 2 + 



1-c 



(4.3) 



(4.4) 



with c = 7V/2 + 1 and z = —s 2 j 2 _ Hn /(AR 2 ). Here I denote the modified Bessel 
function 

I u {z) = e-™! 2 J v {e™l 2 z). 
We obtain 

c NH 2 H + 1 / 2 R H 



bL(s) = 



J 2 -H(h-H,n)j?+H, n nN/2) 

2 (^2)/2 r(Ar/2) 



J(N-2)/2JR 1 jl-H,nS) _ 
(R-ljl-H,nS)( N -W ~ . 



(4.5) 



8 http : //functions . wolfram . com/07 . 34 . 17 . 0012 . 01 
9 http : //functions . wolfram . com/07 . 34 . 16 . 0001 . 01 
10 http : //functions . wolfram. com/07 . 22 . 26 . 0004 . 01 
u http : //functions . wolfram, com/07 . 22 . 03 . 0030 . 01 
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Continue with the remaining case of m > 1. Using (1.4), we obtain 



L( s ) = / c NH s u i G X ' A I - 



N/2 + H 




/2u 



RJu+l(ju,n) 

Using (3.3), rewrite this formula as 



J V {R 1 3v, n u)du. 



,V2i 



2H—m r°° 



RJv+l{jv,n) JO 



L(s) = -5-? — 7i — r / u M R ^nWjGi i 



N/2 + H 




<fu. 



To calculate this integral, use (4.1) with a = to — H + 1, b = R l j v ,n, m = p 
q = l,r = 0,u> = s~ 2 , ai = N/2 + H, and &i = 0. We obtain 



LW = 



CAT// 



nm-ff+1/2 2H-m nra-fl 



1 - 





--H + 1 



X ^3,1 I 2-2 
\° Jv,; 



m-H+l + v N , tt i _ 

2 > 2 """ ' 2 



m-ff+l-f 



To decrease the order of the Meijer G- function from 4 to 2, put v = m—l — H 
and use (4.2). We get 



4R 2 



H+l "2,0 I 2,2 

V * Jm-l-H,', 



1-m + H, N/2 + H 



h i ( \ _ c nh2 t ' s H ^ ,i 

u m — H\Jm.— l — H,n)J rrl —\ — }j n 

For further simplification, use the symmetry relation (4.3). We obtain 



Om-H+l/2 a 2H-mr>m,-H ( q 2 -2 

w /■ \ _ CW-ff^ s £<• ,-,1,0 / * Jm-l-H,; 

"mn\ s ) ~ — — s . m -H+l "0,2 



•Jm-H\Jm-l-H,n)3 m -i-H n 



4R 2 



m-H,l- N/2 - H 



Then, use the argument transformation (4.4) with a = rn/2 + 1/2 — N/4 — H to 
get 



b mn ( s ) 



c NH 2 H +( N -V/ 2 R H +( N ~V/ 2 

S^-^J m -H{j m -l-H,n)C-i% 



xG, 



1,0 / S Jm-1-H,n 
4i? 2 



i/2 + (AT - 2)/4, -m/2 - (TV - 2)/4 



Finally, use (2.4) to obtain 

a , , c NH 2 H +^ 2 R H 



b mn (s) 



J m -H(jm-l-H,n)Cll_ H , n r(N/2) 

v q(;V-2)/2tV ju/^ J m+(N-2)/2{R~ jm-l-H, n s) 



10 
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The last formula and (4.5) can be unified as 

u , s c NH 2 H +^R H 



,(*) 



J\rn-l\-H+l{j\m-l\-H,n)Ji„^l\-H,n T ( N / 2 ) 

x [g m {RT l j\ m -i\-H,ns) - Sq] 



with 



5m (z) = 2(-- 2 )/ 2 r(^/2)^±^f|M. 

Substituting the value of Cnh from (1.2), we obtain 

2 H+1 v / 7r (JV-2)/ 2r ( iv /2 + H)T(H + 1) sm(irH)R H 



L(«) = 



r(JV/2)J| m _ 1 |_ H+1 (j| m _ 1 |_H,n)i| m _i|_ fl -,n 

X [g m {RT l i\ m -i\- H ,ns) - 8^]. 
This result was proved by Malyarenko (2008) for the case of R = 1. 
Example 2. For simplicity, put R = 1. The functions 

¥>„(«) = V2ra+lP n (2u-l), n>l, 

where P n (x) are Legendre polynomials, form a basis in the Hilbert space L 2 [0, 1]. 
The corresponding Fourier coefficients b mn (s) have the following form: 

■ 2 N/2,1 



b l n( s ) = CnhV^u + 1 / II 



H - 1/2 GH I u 



— n, n + 1 
0,0 



G: 



2,0 / u 

2,2 



0,1 -if 



<i-u, 



,(s) = c N hv2ti + ls 



2H-1 



m — tl — 



H-l/2^2.0 
Lx 2,2 



— n, n + 1 
0,0 



: G ia ( ^2 



N/2 + H 




<iu. 



This follows from (1.4) and (2.3) with A = 0. 

The integrals in the last display are complicated, because the arguments of 
the two Meijer G-functions contain different powers of the independent variable 
u. However, they still can be calculated analytically, using formula 2.24.1.1 from 
Prudnikov et al (1990). The answer is 

c NH y/2n + 1 



CM 



(s) 



xG 



<>.(> 



1-2H 3-2H 1-2H 3-2H N 



1 



0,1 -J? 



4 ' 4 ' 4 ' 2 ' 
l-2H+2n 3-2H+2n -l-2H-2n l-2H-2n 

4 j 4 ) 4 ' 4 



c NH \j2n + Is 



2H-m 



x G\i | s 



l-2m+2H 3-2m+2H l-2m+2H 3-2m+2H N_ 

' 2 



4 ' 4 ' 4 

\-2H+2n-2m 3—2H+2n+2m 

4 ' 4 ' 



H 



-l-2H-2n-2m l-2H-2n-2m 



The details are left to the reader. 
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5. Proofs of technical lemmas 

Proof of Lemma 1. It follows from (1.1) that the autocorrelation function of the 
multiparameter fractional Brownian motion can be written as 

R(s, t, u) = l -{s 2H + t 2H - (s 2 - 2stu + t 2 ) H ). 

Therefore, the autocorrelation function of the stochastic process X l m (t) has the 
form 

Rm(s,t) = h - lim h(a), 

a-*(N-l)/2 

where 

c 



h = W S 2H + t 2 H) I C (N-2)/2 {u){1 _ u 2)(AT-3)/a ^ 

j a (a) = C -(2st) H f^u+lT-^l-ur-^ 2 (^±l- u y C ^-^(u)du. 

To calculate Ji, we use formula 2.21.2.17 from Prudnikov et al (1988): 

V - * 2 ) A - 1/2 ^(z/a) dx = flV^ ^ffi 

with a = 1, A = (N — 2)/2, and n = m. Here, ^ is the Kronecker's delta. After 
simplification, we obtain 

, _ 2 N -^( N -^ 2 T((N 2)/2)T((N - l)/2) 2ff ^, 2ffu o 

1 r(jv - 2)r(JV/2) l j m " 

This expression can be further simplified using the doubling formula 12 

T(2z)^ 2 ^T(z)T(z + l/2) 
with z = (N- 2)/2. We get 

TrN/2 

Jl "f(iv72) (s +t )(W 

To calculate ^(a), we use formula 2.21.4.15 from Prudnikov et al (1988): 

(x + a) a_1 (a - a;) A - 1 / 2 (^ - aO _ *C£(sc/a) dx = 

(- 1 )" (Mojl\ n \ OW r(a)r(A+l/2) , > a+A _ 1/2 , , ,_, 
-^(1/2 + A - a)„(2A)„ r(a + A + n+1/2) (2a) (, + a) 

x 3 F 2 (a, 0, 1/2 + a - A; 1/2 + a - A - n, 1/2 + a + A + n; 2a/ {a + z)) 



2 http : //functions . wolf ram. com/06 .02.16. 0006 . 01 
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with a = 1, A = (JV - 2)/2, z = (s 2 + t 2 )/(2st), ?9 = -#, and n = m. After 
simplification, we obtain 

7 2 (a) = (_i)™ 2 tt+3 ( Ar - 3 )/ 2 7r( JV - 2 )/ 2 

r((JV - 2)/2)r((7V - l)/2 - a + m)r(a)r((JV - l)/2) 2H 

r((iV-l)/2-a)r(JV-2)r(a+(iV-l)/2 + m) l J 

/ N - 3 AT - 3 TV - 1 4si 

x 3F2 a, —if, a — ; a m,a -\ \- m; 



2 ' 2 ' 2 (s + «) 5 

Using the doubling formula with z = (N — 2)/2, we get 



\2H 



I 2 (a) = (_i )m2a+ ( J v-3)/ 2 ^( J v-i)/2 r((JV - l)/2 - a + m)r(a)Q, + t)' 
v ; v ; r((AT-l)/2-a)r(a + (iV-l)/2 + m) 

„ / rr iV - 3 N-3 N - 1 4st 

x 3 f 2 [ a, — H , a — ; a m, a -\ \- m; 



2 '" 2 ""'" ' 2 ' ""' (s + t) 2 

(5.1) 
In the case of in = 0, (5.1) simplifies as follows. 

Un\ - o<*+(N-3)/2 (N-p/2 T(a)(s + t) 2H 
21 ■ } T(a+(N-l)/2) 

( N - 3 N-3 N - 1 4s* 

According to paragraph 7.2.3.2 from Prudnikov et al (1990), the value of the gener- 
alised hypergeometric function p F q (a\, . . . ,a p ;b\, . . . ,b q ; z) is independent on the 
order of upper parameters a\, . . . , a p and lower parameters b\, . . . , b q . Moreover, 
formula 7.2.3.7 from Prudnikov et al (1990) states that 

p-T n\(l\ 7 . . . , &p— r: C\ , . . . , C r , 0\ , . . . , 0«_ r , C\ , . . . , C r , Z ) 

- V I -h h ■ \ ^ 

p— r"q— r\U>li • • • i ®"p—ri ^1 •> • • ■ •> ®q— r j % ) ■ 



Using these properties, we get 

>(a) = 2 a +( jv - 3 )/ 2 7 r( Ar - 1 )/ 2 T ( a ^ s + t ^ H F ( c 
21 ' T(a + (N-l)/2) 2 X \ 

In particular, 



) 2H ( N - 1 4st 

i2(«) = 2 7T^ " — ; — 2-^1 ck, — -"! a 



(*+*)* 



lim I 2 (a) = 2 *-V*-i)/ a n(*-l)/2)(« + t) ag 

a->(iV-l)/2 V ^ r(iV-l) 



2 ' ' '(s + t) 2 

The application of the doubling formula with z = (N — l)/2 yields 

lim I 2{a ) = ^-( S + t) 2H 2 F 1 (^-^ 1 -H-,N-l; ^ 



a ^(N-i)/2 zy ' r(A72) v y V 2 ' ' '(* + *) 
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The argument simplification formula 13 states: 

2 F 1 (a, b; 2b; 4z/(z + l) 2 ) = (z + \) 2a 2 F 1 {a, a-b+ 1/2; b + 1/2; z 2 ). 
Use this formula with a = —H, b = (N — l)/2, and z = s/t. We obtain 

n N/2 

lim Ua) = — — r -t 2H 2 F 1 (-H 1 1 - H - N/2; N/2: s 2 /t 2 ). 
a-»(jv-i)/2 2V ' T{N/2) 2 U /»/»/; 

and, finally 

TrJV/2 

i?o(s ' f) = rpj 182 " + <2H(1 " 2jFl(_jff ' : " H N/2; N/2; s2/f2))] - (5 - 3) 

It remains to prove that the first case in (3.2) simplifies to (5.3). To do 
this, use the representation of the Meijer G-function through hypergeometric func- 

-t- 14 

tions 

^■■■' b J h n^=„ + i r(aj - 6 fc ) nj= m+1 r(l - 6, + 6 fc ) z 



x „F„ 



1 - oi + 6fc,...,l - Op + 6/- 



p q 1 yl +oi - a k ,...,l + a k -i -a fe , 1 + a fc+ i - a k) .. .,1 + a g - a fc ; (-l) p 

with m = n = 2, p = g = 4, z = s 2 /t 2 , a\ = 1 — H, a 2 = &i = 0, a 3 = AT/ 2, 04 = 1, 
62 = 1 - H, 63 = 1 - H - N/2, and 64 = — J5T- We obtain 



c 



2 

NH 2H 



1 2 



r(i-ff)r(fl") (h,i,i- N/2,0 



T(N/2)T(N/2 + H)T(1 + H) ^ \^H,H + N/2, 1 + H; s 2 /t 2 



T(H-l)T(2-H) , s 2(i-g) (1,2- H,2-H- N/2,1 -H; 

+ T(N/2 + H- l)T(H)T(N/2 + 1) ^(l-ff) 4 3 ^2 - iJ, JV/2 + 1, 2; s 2 /i 2 

The first term is simplified, using the following formula 15 
p F q (0,a 2 , ■■■ ,a p ;bi,...,b q ;z) = 1, 
while the second term is simplified by (5.2). We get 

j = iff s 2H [ T(1-H)T(H) 

1 2 [T(N/2)T(N/2 + H)T(1 + H) 

T(H-l)T(2-H) gaO-H) f 1,2 - H - N/2,1 - H; 

+ T(N/2 + H - l)T(H)T{N/2 + 1) £2(i-ff) 3 2 ^ iV/2 + 1,2; s 2 /t 2 

Using the formula 16 

3 F 2 (l,b,c;2,e;z) = —-^——[ 2 F 1 {b-l,c-l;e-l;z)-l}. 
{b- l)[c- \)z 

13 http : //functions . wolfram. com/07 . 23 . 16 . 0005 . 01 
14 http : //functions . wolfram. com/07 . 34 . 26 . 0004 . 01 
15 http : //functions . wolfram. com/07 .31.03. 0012 . 01 
16 http : //functions . wolfram. com/07 . 27 . 03 . 0120 . 01 
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with b = 2 - H - N/2, c=l-H,e = N/2 + 1, and z = s 2 /t 2 yields 



c: 



2 

NH „2ff 



1 2 



r(l - H)T(H) 

T(N/2)T(N/2 + H)T(1 + H) 

T(H - 1)T(2 - H)(N/2) t 2H 

+ T(N/2 + H - l)T(H)T(N/2 + 1)(1 - H - N/2)(-H) s^? 
x [ 2 Fi(l - if - AT/2, -H; N/2; s 2 /t 2 ) - 1]] . 

To simplify the second line, use (2.1) in the following way: 

r(H-l) = T(H)/((H-l), 
T(2-H) = (1 -H)T(1 -H), 
T(N/2 + H - 1)(1 - H - N/2) = -T(N/2 + H), 
T(H)(-H) = -T(H + 1), 

T(N/2 + 1) = (N/2)T(N/2). 

After simplification, we obtain 

'■ = ^/trwn™ + 1) [""-""(i-^ic-g. x-H-m m w>\. 

or, by (5.4) with z = H, 
h 



c NH n 



2T{N/2)T{N/2 + H)T(H + 1) sm(irH) 

x [s 2H + t 2H (l - 2 F 1 (-H, l-H- N/2; N/2; s 2 /t 2 ))}. 

Substituting (1.2) to the last display, we get 

TT N / 2 

I i = f ^j Y) [s 2H + t 2H (l- 2 F 1 (-H,l-H-N/2;N/2;s 2 /t 2 ))}. 

This completes the proof of the case of m = 0. 

In the case of in > 1, (5.1) can be rewritten as 

lim /,(«) = (-i r2 "-y*-i>/» (m - 1)W ^ ~ 1)/2 \ (g + ^ 
«-(jv-i)/2 v ' y ' Y{N-l + m) 

3 F 2 (a,-H,a-^;a-^-m,N-l + m; ^p) 

x lim 



a-(Ar-i)/2 r(a - (JV - 3)/2 - m) 

r(a-(iV-3)/2-m) 

x lim — — r-^ -r 1 r — - . 

a ^( N -i)/2 T{{N - l)/2 - a) 
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To calculate the first limit, we use the following formula 17 

p Fg(a 1 ,...,a p ;bi,. ..,b q ;z) _ z n+1 IIf=i( a j)n+i 
bl ™« T(h) ~(»+l)!ni=a(6iWi 

x pF q (ai + n+ 1, . . . ,a p + n+l;b 2 + n+l, . . . ,b q + n+ l } n + 2;z), n G Z + 

with n = m — 1, p = 3, q = 2, ai = (JV — l)/2, a 2 = — -ff, a 3 = 1, b\ = 1 — to, 

b 2 = N - 1 + to, and z = 4st/(s + i) 2 . We get 

JT f™ W ™ W-3 . „, iV-3 

5^2 a, —-ti,tt — 



3 r 2 I a, — -n,a g - ' u 2 m > iv ~~ i "T "'» (a-ft) s 

«_,( ™i)/2 r(a - (JV - 3)/2 - to) 

2 3m (st)"T((iV - l)/2 + m)r(m - H)T(N - 1 + to) 

(s + t) 2m r((7v - i)/2)r(-i?)r(iv -1 + 2™) 

x 3 F 2 ((iV- l)/2 + TO,TO-.ff,TO + l;iV- 1 + 2to,to + l;4si/(s + £) 2 ). 
Using the doubling formula with z = (N — l)/2 + to and (5.2), we get 

3 F 2 (a, -H, a - ^; a - ^ - to, JV - 1 + to; ^ 

lim 

a-»(w-i)/2 r(o-(JV-3)/2-m) 

(st) m T(m - H)T(N - 1 + m)y/7T 

(s + t) 2m T((N - l)/2)T(-H)2 N - 2 T(N/2 + m) 

x 2Fi((N - l)/2 + to, to - #; TV - 1 + 2to; 4st/(s + t) 2 ). 

Rewrite the second limit as 

r(q - (JV - 3)/2 - to) _ r(ot - (TV - 3)/2 - m) 

a-(iv™i)/2 r((iV- l)/2-a) ~ a-( ™i)/2 r(a - (iV- l)/2 

x i im r(a-(iV-l)/2) 
a-(JV-l)/2 r((iV - l)/2 - a) " 
For the first part, use the formula 18 



T(z-n) 



n 



Viz) 1J -2-Jfc 

v ; fc=i 

with z = a — (N — l)/2 and n = to — 1. For the second part, use (2.1) with 

z = a - (TV - 3)/2. We get 

r(a - (TV - 3)/2 - to) = (-I)"' 1 (-1)™ 



o-(jv-i)/2 r((7V - 1)/2 - a) (to-1)! v ' (m-1)! 

Combining everything together, we obtain 

7 r Ar / 2 (s + f) 2 ( iT - m )(sf) m r(m - H) 



Rm(s,t) 



T(-H)T(N/2 + m) 
x 2-Fi((7V - l)/2 + m,m-H;N-l + 2m; 4st/(s + t) 2 ). 



17 http : //functions . wolf ram. com/07 . 31 . 25 . 0003 . 01 
18 http : //functions . wolfram. com/06 . 05 . 16 . 0022 . 01 
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Use (2.1) with z = 1 - H: 

tt n / 2 (s + t) 2 ( H - m \st) m T(m - H)H 
m(M) ~ T(l-H)T(N/2 + m) 

x 2 ^i((A r - l)/2 + to, m - H; A - 1 + 2to; 4s*/ (s + t) 2 ). 
By the additional formula for the gamma function 19 

r(z)r(i - z) = -^— 

sm(7rz) 
with z = H, and (2.1) with z = H + 1 we get 

n (N-2)/2( s + t)^ H - m ){ s t) m T{m - H)T(H + 1) sin(Trff) 

T(N/2 + to) 

x 2 ^i((A r - l)/2 + to, to - H; A - 1 + 2to; 4s*/ (s + i) 2 ), 
and then, using the following formula 20 

(V5 + .)-/,W; ! M«(Vl + 1) 2 ) = r(> + i/2)r«,- a + i/2) 

r(o) 

l-o,6-o + l/2 

0,1/2-6 



(5.4) 



R m {s,t) = 



X 6*2^2 I : 



with z = s /t , a = to — H, and b = (A — l)/2 + to, we finally obtain 



R m (s,t) 



NH „m.2if-m/il,l 



S"t 



G " • ' *^ 



r 2,2 



if + 1 - TO, A/2 + if 

0, 1 - AT/2 - m 



m > 1. 



□ 



G 2 ; 2 I ■'• 



This completes the proof. 

Proof of Lemma 2. (3.3) for the case of to = is obvious from the formula 21 

r(oi + o 2 - 6i - b 2 ) 

x 2-Pi(o2 ~ 6 2 , Oi - 6 2 ; oi + a 2 - &i - 6 2 ; 1 - a;) 

with x = u 2 /s 2 , ai = iV/2, a 2 = 1, &i = 0, and b 2 = 1—H, where we use a shortcut 
(1 — x)+ for max{l — x, 0}. Similarly, (3.3) for the case of to > 1 is obvious from 

X b (l - x) a + b ~ l 



the formula 



22 



-,1,0 



G{' U i * 



T(a-b) 



with x = u 2 /s 2 , a = A/2 + H, and 6 = 0. 

It remains to calculate two integrals. The first one is as follows: 



h 



-NH 



i{».t} 



u 2H - l G 2 >° 



o 



2,2 I 9 



A/2,1 
0,1 -H 



G 



2.2 



f 2 



A/2,1 
0,1- J? 



du. 



19 http : //functions . wolfram. com/06 . 05 . 16 . 0010 . 01 
20 http : //functions . wolfram. com/07 . 03 . 26 . 0031 . 01 
21 http : //functions . wolfram. com/07 . 34 . 03 . 0645 . 01 
22 http : //functions . wolfram. com/07 . 34 . 03 . 0247 . 01 
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Taking into account (3.3), we can substitute oo for the upper limit of integration. 
After change of variable u — s/x, we obtain 



h 



~NH 



H-l r 2,0 
J, <~ r 2 2 



il 



N/2,1 



Gi 



o,i-ff 2 ' 2 * 2 



N/2,l 
0,1 -H 



dx. 



To calculate this integral, use (2.2) with a = H,s = m = u = v=p 

t = v = 0, w = 1/s 2 , z = 1/t 2 , ci = Oi = N/2, c 2 = a 2 = 1, d\ = b\ - 
d 2 = b 2 = 1 - H . We get 

, = ifl2M 2 ^ 1-H,0, N/2,1 
1 2 4 > 4 \t 2 0, 1 - J?, 1 - H - N/2, -H 

This completes the calculation of the first integral. 
The second integral is as follows. 



= 9 = 2, 
0, and 



-NH 



(St) 



2H-m 



a{s,t} 



n ,2m-2H-l/^l,0 

^1,1 \ s 2 



N/2 + H\ n i,ofu 







G M 



Vt 2 



N/2 + H 




Taking into account (3.3), we can substitute oo for the upper limit of integration. 
After change of variable u — -,/x, we obtain 



h = 



^NH 



(St) 



2H-r 



i-H-l n lfi I x 



N/2 + H 







1,0 / ^_ 
t 2 



N/2 + H 




dx. 



To calculate this integral, use (2.2) with a = m — H, s = m = u = v = p=q=l, 
t = v = 0,w= 1/s 2 , z = 1/i 2 , ci = oi = A^/2 + if, and d x = 6 X = 0. We get 



-NH m+2H-mfil,l 



2 H+l-m, N/2 + H 

0, 1 - N/2 - H 



This completes the proof. 



□ 



du. 
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